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Abstract 

Q Upward throughgoing muons in underground detectors are induced 
by atmospheric with energies in the 1 — 10^ GeV range. The oscil- 
lations of atmospheric muon neutrinos should affect mainly those with 
energy smaller than E^, ~ 50 GeV, or < 10 GeV . We analyzed 
the MACRO0 upward through going muons data set, using the Multiple 
Coulomb Scattering to estimate the energy E"^ of muons crossing the 
detector. We analyzed the event distribution for six event subsamples 
with different values of < ly/Ey >. The event distribution is in agree- 
ment with the hypothesis of neutrino oscillations with Am^ few times 
10"'^ eV'^ and maximum mixing. 



1 Introduction 

Neutrino oscillations, in a two neutrino mixing scenario, are the most likely 
solution for the atmospheric neutrino problem Underground experiments 
unfold the mass difference Am^ = m^j^g — m^j,^ and the mixing angle 6 from 
the measurement of the survival probability 

The survival probability is usually plotted versus the measurable quantities ii, 
and/or E^,. i^, is the neutrino path length, estimated from the measurement of 
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the zenith angle B of the daughter lepton. The neutrino energy E^, is estimated 
from the calorimetric measurement of the energy of all final state particles. 

At present, E^y is measured with some accuracy for low energy events (fully 
and partially contained SK [Q] and Soudan2 P] events); however in this case ii, 
is poorly reconstructed. The vice versa happens, when in the final state only 
the high-energy muon is measured (high energy neutrinos interacting in the 
rock around the detector, and inducing up-throughgoing muons in MACRO 
[0], SK and Baksan^]). The goal of next generation of atmospheric z/ exper- 
iments should be the measurement of both and E,y, with enough precision 
to observe the first minimum of ^(t'^ v^x)- 

In the case of the MACRO experiment, the measured oscillations parame- 
ters (best value: Am? = 2.5 x 10~^ eV^ and sin'^26 = 1) are deduced: from the 
deficit and the distortion of the zenith angle distribution of up throughgoing 
/i {i.e. -P(t'^ — > ty^) vs. cosQ) 0]; from the measured deficit of lower energy 
partially contained events^]. 

In this paper we present the following approach: we identified topologies 
of upgoing muons induced from neutrinos with different energies and path 
length and we discuss the results in terms of oscillations. The energy 
E^ of detected muons is estimated through the multiple Coulomb scattering; 
Efx is correlated with the energy E^, of the parent neutrino. The survival 
probability shows a dependence with the < ly/Ey > parameter, unexpected 
without neutrino oscillations, while it is in good agreement with the quoted 
oscillation parameters. 

2 Up throughgoing muon energy separation 

Up throughgoing muons are induced both by (relatively) low energy atmo- 
spheric z/^, interacting near the detector, and higher energy neutrinos, inter- 
acting farther. Assuming the neutrino oscillations mechanism, with Am^ = 
2.5 X 10~^ eV"^ and sin'^26 = 1, a reduced flux of muons with E^ < 10 GeV is 
expected, as shown by the shadowed histogram of Fig. la. 

Multiple Coulomb Scattering. Muons traversing the detector are deflected 
by many small-angle scatters, the bulk of which are Coulomb scatterings from 
the nuclei and from the atomic electrons (Multiple Coulomb Scattering, MCS). 
As an effect, muons are slightly deviated from a straight-line direction by an 
angular deflection (projected in a plane) (ppiane 0- Quantitatively, the averaged 
angular deflection (pp^^e spatial deflection y^'^^ in one projective plane 

are approximated with: 
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Figure 1: a) upper: energy spectrum of simulated neutrino- induced up throughgoing 
in MACRO. Full histogram for no oscillations; shadowed histogram, for {v^ — > v^) 
oscillations, with Am^ = 2.5 x 10^^ eV'^ and si'n?29 = 1. b) lower: profile of the 
averaged value of the parameter Vyj vs. muon energy. The full (open) points refers 
to events generated with the full (table) MC simulation, see text. 

In the case of MACRO, whose radiation length is Xq = 22.6 gr/cm^ and 
total thickness x = 400 gr/cw?, the average angular deflection is (jf^ne — 
0.052 {rad)/E^ [GeV). This corresponds to a projected spatial deflection in 
the wires view of the limited streamer tubes (whose cell size is 3 cm) of y'^'^'^ ~ 
7{cm)/E^{GeV). 

Experimental technique. The technique we used to distinguish muons of 
low and high energy is illustrated in Fig. 2. All up throughgoing /i, crossing the 
whole apparatus, were analyzed. The streamer tube hits in the lowest planes 
(N=l:8) were used for a refit, and the straight-line fit parameters were used to 
define the lower track. In a similar way, the highest planes (N=6:14) were used 
to define the upper track. The upper and lower tracks are used to evaluate r^, 
the spatial difference in the 2 = plane between the two tracks, and A$, the 
angular difference between the slopes. The value of the parameter {^^) 
depends on the muon energy E^. In Fig. lb the profile [i.e. the average value 
of evaluated in each bin of E^ from our MC events) is shown. From the 
figure, the spatial difference between the lower and upper track exceed the 
intrinsic resolution of streamer tube cells when E^ < 10 GeV. The behavior 
of the A$ parameter is analogous. 

Event classification. A selection cut, based on r^^ and A$, was defined. Each 
up throughgoing muon is classified as: L= Low energy if > 3 cm and A$ > 
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Figure 2: Sketch of the refit procedure to obtain a lower and an upper track, and 
the definition of the parameter r^ . 

0.3°; H= High energy if < 3 cm and A$ < 0.3°. The remaining events with 
ambiguous deviation, were classified as M=Medium energy. Each sub sample 
has a different average energy for the parent neutrino. If z/^ disappearance is 
at work, the reduction should affect mainly the L sample, and we expect in the 
data a reduced ^ ratio, with respect to the MC prediction without oscillations. 

Checks on Monte Carlo predictions. A crucial point is the MC reliability, 
to avoid any systematic bias which can invalidates the results. For a detailed 
check of the detector simulation, we compared two large samples of measured 
with simulated downward going atmospheric muons. Down throughgoing at- 
mospheric muons have an average residual energy of ~ 300 GeV when reaching 
the underground laboratory, and 90% of them have > 10 GeV. Our cut 
classified 63% (25%) of these events as High (Low) energy. The small frac- 
tion of atmospheric muons (~ 0.3%), reaching MACRO with ii^^ < 1 — 5 GeV, 
stops inside the detector: 88% (7%) of them have been classified as Low (High) 
energy. 

The energy distribution of these two atmospheric samples overlaps the spec- 
trum of neutrino-induced upgoing /i: the down-stop events in the lower, and 
the down- throughgoing in the high-energy tail. Using atmospheric muon sam- 
ples, we adjusted the streamer tubes geometry in the MC database. After the 
MC correction, the measured and expected distribution of and A$ variables 
match together, both for the down-stop and down-throughgoing atmospheric 
muons. 
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3 Results 



Real data. We used the same up throughgoing muon sample described in 
[|l|]a. The details of the events identification, cuts and selection efficiency were 
published 0, In our analysis, only the 316 events crossing the whole appa- 
ratus where used, to enhance the MCS effect on the muon track. The and 
A$ parameters are evaluated from the reconstructed lower and upper tracks, 
and 101 events were classified by the cuts as Low, 51 as Medium and 164 as 
High energy. 

Simulated data. The expected number of up throughgoing muons from 
interactions was evaluated with two different methods. In the first method 
("table" MC) 0, ^, the muon energy spectrum is obtained from a numerical 
convolution of the differential atmospheric neutrino spectrum, neutrino cross- 
sections and muon propagation in the rock. In the second ("full" MC), a 
full simulation is used: the z/^ interact in a large volume of rock around the 
detector, and the muon is propagated through the rock up to the detector. 
The simulated data have the same format of the real ones, and thus followed 
the same analyses steps. We applied the cut on and A$ on the selected up 
throughgoing MC muons: for a live time period equivalent of that of the real 
data, we expect 431 (430) events from "table" ("full") MC: 177 (178) Z, and 
191 (193) H. 

Uncertainties. The overall theoretical uncertainty on MC predictions is 17% 
[1,2]. This is a scale factor, which is almost constant when the number of events 
is plotted versus cosG or E^. The shape of the distributions (as the angular 
one presented in @, 0) has a reduced theoretical uncertainty: we assumed 
a (conservative) value of 5% in the considered range of v energies. The sys- 
tematic uncertainties depend from many factor: the number of upper/bottom 
planes used in the refit; the choice of a particular value of and A$; the 
MC simulation statistics; the fluctuation in the streamer tube and scintillator 
efficiency along the running period; the detector acceptance uncertainty vs. 
the muon zenith angle. The evaluated overall systematic uncertainty is 7%. 

Data and MC comparisons. We divided the events in three bins of 
{L,M,H). The L and H samples were sub-divided in two bins of neutrino 
path length variable ii, (the Medium energy sample have too few events). The 
number of detected and expected events vs. < logiQ^i^/ E^) > (MC evaluated) 
is presented in Table 1 and Fig. 3. It is also included an additional point, 
corresponding to the partially contained lU events 0. These events are in- 
duced by neutrinos with an average energy E'^, ~ 4 GeV, interacting inside the 
detector. This sample is contaminated by a 13% contribution of z/e and NC 
interactions. 

The ratio R between the detected to the expected number of events for 
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Figure 3: Ratio R of tlie measured to the expected number of events vs. < 
logio{i^/E^) >. The H and L samples are subdivided in two intervals of i,^. The 
black triangle is the point corresponding to the partially contained upgoing events 
is included [4]. The horizontal full line represents the expectation assuming no 
oscillations. The segment below each data point is the expectation assuming — > 
i^x, with oscillation parameters quoted in Fig. la). 
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Figure 4: a) Full histogram: simulated distribution of reconstructed up through- 
going muons in MACRO vs. logio{£u/E^). The shadowed region represents the ex- 
pected distribution in case of neutrino oscillations (with the parameters of Fig. la), 
b) The full histogram corresponds to the distribution of simulated neutrino-induced 
muons, classified of high-energy from the Multiple Coulomb Scattering analysis, and 
with large zenith angle. The dashed histogram is for low-energy muons, with small 
zenith angle. The average value of these two distributions, are plotted in Fig. 3 
(points 1 and 5). 
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Table 1: Number of detected and expected events {N(iata= detected, Nmc= ex- 
pected and Ar^g=expected with oscillations). In column 3, there is the average 
value ± the HWHM of the logiQ{£i, / E^) distribution. 



each sub sample is shown in Fig. 3. The value of R expected with the used 
flux and neutrino cross sections, is 1 (no oscillations). This value has an overall 
scale factor uncertainty of 17%, which is almost constant all over the {li,/Ei,) 
interval. In our data, R decreases when < logiQ{£i, / E^,) > increases. This 
behavior is consistent with the neutrino disappearance probability, when folded 
with our distribution of events. In Fig. 3 the segment below each data point is 
obtained from Monte Carlo simulated events, assuming oscillation parameters 
Am^ — 2.5 X 10~^ eV^, sin^29 — 1. The effect of neutrino oscillations in the 
distribution of {Cy/ Ey) is shown in Fig. 4a as a shadowed histogram. The effect 
of our cuts on the simulated events is shown in Fig. 4b for two subsamples 
(the one with the lowest /highest average value of {iy/E^)). 

We evaluated the foi" the shape of the ratio R vs. < logio{iiy/Ey) > for 
the two different test hypotheses: T0= no and Tl= neutrino oscillations with 
the above quoted parameters. For hypothesis TO, we expect a (almost) con- 
stant value. Normalizing the number of expected with the number of detected 
events (316/430), we get a x^/dof = 14.8/5 (probability of less than 2%). For 
hypothesis Tl we expect a value of R which decreases as a "ladder" from 1 
to 1/2, when {e^/Ey) > 1000 Km/GeV. In this case, we get x^/dof = 4.7/5 
(probability=45%) 

Hypotheses TO and Tl can be tested also using the ratio of ratios {jj)data/ (^)mc- 
With the double ratio, most of the theoretical and systematic uncertainties 
cancel. For Tl, we expect a larger reduction of events with Low energy with 
respect to the ones with iJigh energy. From our data, {jj)^^^ = 0.62 ± 0.08. 
For TO, {^YiSc = 0.92 ± 0.08 (c.l. of 0.5%). For Tl {jjf^c ="0.68 ± 0.06 (c.L 
90%). 
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4 Conclusions 



We discussed the possibility to study the present high-energy atmospheric 
neutrinos data in terms of different bins of < logioi^ly / Ey) > , using the 
Multiple Coulomb Scattering on muons. The MACRO events were analyzed: 
the result is consistent with the hypothesis of neutrino disappearance with 
Am^ = 2.5 X 10~^ eV"^ and sin^29 = 1. This strengthens the hypothesis 
of neutrino oscillations to explain the data, against alternative more exotic 
models. A global analysis of all MACRO neutrino-induced data sets (which 
includes the results of this work) is in progress. 
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